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Abstract 

By using a hydrodynamic atmospheric escape mechanism (Levi and Podolak 



2009 1 we show how the unusually high mass density of Quaoar could have been 
predicted (constrained), without any knowledge of a binary companion. We 
suggest an explanation of the recent spectroscopic observations of Orcus and 



Charon (Delsanti et al. 2010 Cook et al. 2007). We present a simple relation 



between the detection of certain volatile ices and the body mass density and 
diameter. As a test case we implement the relations on the KBO 2003 AZg4 
and give constraints on its mass density. We also present a method of relating 
the latitude-dependence of hydrodynamic gas escape to the internal structure 
of a rapidly rotating body and apply it to Haumea. 

Key Words: ICES; KUIPER BELT; TRANS-NEPTUNIAN OBJECTS 
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1 Introduction 

Jeans escape is often used to estimate the escape rate of volatiles from KBOs (see 



Schaller and Brown 2007| for example). However, for small bodies hydrody- 



namic escape may be more important. A good discussion of atmospheric escape 



is given in Stern and Trafton ( 2008 1 . In previous work ( Levi and Podolak 2009 



hereafter paper I) we have shown that, with proper modifications, the coronal 



flow theory, as derived by Parker (1963), can be used to describe hydrodynamic 



atmospheric escape off some KBOs. We found that for an adiabatic hydrody- 
namic atmospheric escape to occur, the body's diameter D has to be smaller 
than some maximum value, D maxi which obeys the relation: 



y TT&p b m g 7-I 

Here k and G are Boltzmann's and Newton's constants respectively, Tf, is the 
body surface temperature, p b the body mass density, m g the molecular mass of 
the atmospheric component of interest, and 7 its adiabatic index. 

It is important to note that the mechanism of volatile loss and retention for 
KBOs in hydrodynamic flow is inherently different from that of Jeans escape 



( Schaller and Brown 2007 1 . An in-depth discussion of the differences between 
the two mechanisms can be found in paper I. We simply point out here, that, 
although the Jeans mechanism may indeed control atmospheric escape from 
certain large KBOs, it is often a slow process. Although a low rate, integrated 
over the age of the solar system can remove the required amount of volatiles, a 
rapid replenishment by a short-term process such as cryovolcanism will deposit 



volatiles on at least part of the surface before Jeans escape can fully remove 
them. The assumption of Jeans escape must be coupled with the assumption 
of no fast replenishment otherwise it can lead to erroneous interpretations. A 
hydrodynamic escape mechanism removes gas much more quickly, and so avoids 
this limitation. A case in point is Enceladus, where jets of water emanate and 



produce a plume of gas containing CO and other gas species ( Waite et al. 2009 1 . 



However CO ice is not found as a surface deposit near the south pole (Brown 



et al. 2006 ) because it escapes before it can recondense. 



It is commonly assumed that unprocessed water ice in outer solar system 



objects is amorphous (Notesco and Bar-Nun 1996 Prialnik et al. 2008 Je- 



witt et al. 2008). However, crystalline water ice was observed on the surfaces 



of many intermediate size KBOs, e.g. Charon (Cook et al. 2007), Haumca, 



2003 AZ 8 4, Orcus, and Quaoar [see Delsanti et al. (2010) and references cited 
therein]. The presence of crystalline rather than amorphous ice suggests that 
the surface must have been heated to a temperature of approximately 100 K, 
far above the 30 — 40 K surface equilibrium temperatures typical of KBOs. The 
surface temperature of Pluto, for instance, has been put at about 38 K based on 



surface vapor pressure estimates (Delsanti and Jewitt 2006). Although Porter 



et al. (2010) have suggested that micrometeorite impact annealing may convert 



amorphous to crystalline ice in the outer solar system, it is not yet clear if the 
micrometeorite flux in the Kuiper belt is sufficiently high to process the ice 
on those bodies. In any case, the precise heating mechanism is not important. 
What is important is that the surface was heated to at least 100 K for some 
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period of time. In this paper we consider how these results may be used to 
estimate the density of KBOs. In section 2 we present the basic assumptions 
underlying the model. In section 3 we consider the relevant timescales. In sec- 
tion 4 we apply this model to several KBO's for which the required observational 
data are available. Here we also extend the model to rapidly rotating bodies. 
In section 5 we discuss the results. 

2 Basic Considerations 

Certain volatile species that are mostly frozen at 30 — 40 K, will experience 
fast desorption when the temperature reaches w 100 K, establishing a transient 
gaseous atmosphere. If the diameter of a body is larger than D max it will have 
a strong enough gravitational field to prevent the gas from flowing away hydro- 
dynamically. The gas will then either recondense upon re-cooling of the surface 
to the equilibrium temperature or flow via tangential winds to colder regions 
where, again, recondensation would ensue. On the other hand, if the diame- 
ter is smaller than D max the volatile atmosphere will escape hydrodynamically. 
In this case repeated heating of different surface ice patches may eventually 
clear the surface of the volatiles considered. If the body allows a particular gas 
species to flow hydrodynamically into space then that species will not appear 
in the body surface spectra. We now ask for which volatiles such a mechanism 
will apply. Note that wc do not require the loss of large amounts of volatiles. 
All we require is that if surface volatiles are vaporized, they be able to leave the 



G 



body before they recondense on its surface. 

For convenience, in what follows, we refer to cryovolcanism as the heating 
mechanism. However, we have made no effort to identify the precise mecha- 
nism or to estimate its strength. The fact that the ice was able to undergo the 
transformation from amorphous to crystalline form, means that some mecha- 



nism was able to provide the necessary energy. Studies by Schmitt et at. ( 1989 \ 



have shown that the time for amorphous ice to transform is very temperature 
dependent. At T — 100 K it takes hundreds of years, while at T = 140 K it 
takes on the order of an hour. If the methane is in the form of a fairly homoge- 



neous layer, then its vapor pressure will be ~ 4 x 10 Pa at T = 100 K (Brown 



and Ziegler 1980). A surface ice with a vapor pressure given by P vap (T) will 



sublimate with a flux of (Levi and Podolak 2009) 



T = P v , 



2nkT 



9 "cm'V 1 



If pi ce is the bulk density of the ice, then the thickness, £, of the ice layer changes 
at a rate 

d£ _ T_ 
dt p 

The density of solid methane at cryogenic temperatures is of the order of 
0.5gcm~ 3 (Bol'shutkin et al. 1971), and this gives a vaporization flux of 7 gcm~ 2 s _1 . 



If the surface remained at this temperature for the 100 yr required for amorphous 
ice to transform to crystalline ice, a layer more than 4 x 10 5 km thick would be 
lost. At T = 140 K, the vapor pressure is 1 x 10 6 Pa, and the vaporization flux 
is some 160gcm _2 s _1 . In the hour it takes the amorphous ice to transform, 
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around 11 km of ice would be lost. Any water mixed in with the methane would 
remain frozen and unless the methane layer was very pure, an ice layer would 
quickly be formed that would cover any remaining methane. If the methane is 
adsorbed on a water ice layer initially, the vapor pressure above the ice would 
be considerably lower. We consider this case below. 

In what follows we assume that the volatile in question is adsorbed on the 
water ice. In general, the bond between a volatile and water ice is stronger than 
between the volatile and itself, so this assumption gives an upper limit on the 
escape time. We further assume that the dcsorption timcscale is the limiting 
factor for escape. Because we are considering hydrodynamic flow, once the gas 
enters the atmosphere and is allowed to flow hydrodynamically, it will escape 
the body in a very short time (of the order of hours) . We also assume that areas 
where crystalline ice is observed were heated to at least 100 K. If the body is 
small enough, the vapor will flow adiabatically, as described in paper I, and it 
will escape hydrodynamically. The flow will be essentially radial and the gas will 
not migrate to some colder (shielded) region. Finally, we argue, on the basis of 
symmetry, that we expect the KBO surface to be heated more or less uniformly. 
It therefore seems likely that if crystalline ice is observed on the surface of such 
a body, we should not expect large shielded areas where surface volatiles have 
not been exposed to heating. However, if such shielded areas existed, and were 
observed, and if the amorphous ice in these areas were covered by methane or 
some other material, then our argument would not apply. 
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3 Timescales 



Suppose that volatiles are replenished by short-term processes such as cryovol- 
canism. If the adsorption timescale is longer than the replenishment timescale, 
the volatiles cannot escape between replenishments, and volatile loss by the 
Jeans mechanism is unlikely. Here we estimate the adsorption timescale. The 
average time, r, that a molecule stays adsorbed on an icy surface is given by: 



v N e 



(2) 



Where vm is the frequency of oscillation of the adsorbed molecules and E a is 



Hobbs 


1974 


). As in 


Leger 


(1983) 



atom surface vibrational frequency which is of the order 10 Hz. 

CH 4 has neither a dipole nor a quadrupole moment, thus it is not expected 
to interact strongly with OH bonds at the water ice surface. Adsorption is very 
sensitive to the surface coverage percentage and to the number of favorable ad- 
sorbing locations on the icy surface such as defects and amorphous locations. 
Adsorption measurements of CH4 on water ice I/j result in heats of adsorption, 
E a /k, in the range 1380 — 2050 K, where the lower bound is for ice heated to 



238 K and then re-cooled to 77 K (Chaix and Domine 19971. The heating 



probably causes a better ordering of the ice surface and eliminates defects and 
amorphous structures. The surfaces of KBOs continuously suffer impacts of 



cosmic rays and solar UV and visible photons (Cooper et al. 2003 Cook et al. 



20071 that break water bonds and reduce crystal water ice to amorphous ice. 



Indeed, observations of crystalline water ice on KBO surfaces are usually ac- 



companied by the detection of some amount of amorphous ice as well. We shall 
therefore use the upper bound for the heat of adsorption of methane to water 
ice. 

We found no data for the heat of adsorption of acetylene (C2H2), ethylene 
(C2H4), and ethane (C2H6) on water ice in the literature. However, the heats 
of sublimation, AH su b/k, for the pure volatiles for these hydrocarbons are: 



2827, 2202 and 2466 K, respectively ( |Stephenson and Malanowski| |1987 ). The 
bonding energy of volatiles with water ice is generally stronger than the bonding 
with their own kind, so the heats of adsorption to water ice will most likely be 
higher than the values cited. We therefore feel confident that the values for 
Ejk for CH 4 , C 2 H 2 , C 2 H 4 , and C 2 H 6 fall between 2000 - 3500 K. They are 
unlikely to be much higher than this since no hydrogen bonding exists between 
these hydrocarbons and water. For this range of the heat of adsorption to water 
ice, and for Tb = 35 K we find r ~ 10 3 — 10 22 yr. Note that the upper bound 
greatly exceeds the age of the solar system. 



Cooper el al. (2003) have shown that the timescale of amorphization of 



crystalline water ice, along the depth probed in surface spectroscopy, by galactic 



cosmic rays is of the order of 10 6 — 10 7 yr. Cook et al. (20071 have calculated that 
the action of solar photons in the range 100 — 620 nm is even more effective, with 
an amorphization timescale of 10 4 — 10 5 yr, suggesting a periodic replenishment 
of crystalline water ice on the surface of Charon, through cryogenic activity, on 
time scales of less than 10 yr. Therefore, even our lower bound on r is only 
an order of magnitude shorter than the expected time between cryovolcanic 
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eruptions. If some area on the object's surface, perhaps near the poles, reaches 
temperatures as low as 30 K then the lower bound for r becomes 10 7 yr. Thus 
the hydrocarbons we have mentioned above may be considered frozen on the 
surface (or at least some part of it) between successive cryovolcanic eruptions. 
When a surface patch of ice is heated to w 100 K, the average time of adsorption 
reduces to 10 -5 — 10 s. In such a case the hydrocarbons considered are promptly 
released to take an active part in forming a transient atmosphere. The question 
is whether the body will allow the transient atmosphere to flow away during the 
excess heating period or not. A clue to the answer is contained in eq.[l] 

Methanol is particularly interesting in this respect. Molecular dynamics cal- 
culations give a value of 5992 K for the heat of adsorption of methanol (CH 3 OH) 
on water ice. This high value is due to two hydrogen bonds formed between the 



two molecular species ( Jedlovszky et al. . 2006 1 . Thus even for a temperature of 
100 K the adsorption time would be r ~ 10 4 yr. Since the re-cooling of cryo-lava 



takes less than 1 yr Lorenz (1996), the surface will re-cool before methanol can 
take part in a transient atmosphere. For an ice surface temperature of 35 K 
methanol will have an adsorption time of 10 53 yr, much longer than the age of 
the solar system. So methanol should play no part in a transient atmosphere 
and the body size (i.e. gravity) should play no role in deciding whether methanol 
shows up in the surface spectroscopy. It should only be a question of whether 
some chemical process was able to transform methanol to some other species. 

A tentative corroboration of this argument comes from the possible obser- 
vation of methanol on the small object Bienor which also, apparently, shows 
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crystalline water ice in its spectrum [see ( Guilbert et al. , 2009 ) and references 
cited therein]. This implies that the surface was temporarily heated to about 
100 K. Because of its small diameter of 200 km, Bienor would require a mass 
density greater than 41gcm~ 3 in order to retain methanol. 

CO (and certainly also N2) is highly volatile, having a heat of adsorption 



to water ice of 1150 K (Nair and Adamson 19701. Therefore at surface tern 



peratures of 30 — 40 K CO will have adsorption times of t ~ 10 2 — 10 -2 s. So 
even at KBO equilibrium temperatures, in between the periods of excess sur- 
face heating, CO (and N2) would be a major atmospheric component that can 
either escape hydrodynamically or via a fast Jeans mechanism (paper I) and 
may become depleted in between cryovolcanic eruptions. CO and N 2 therefore 
demand a more rigorous treatment such as that given in paper I, which allows 
for hydrodynamic escape. 

4 Results 

For each volatile of interest, we need to determine the appropriate values of its 
polytropic index, 7. At a temperature of 100 K the translational and rotational 
degrees of freedom are active and behave classically whereas molecular vibra- 



tional degrees of freedom are largely frozen (Loeb 1961). Using experimental 



results for the isobaric heat capacity, C p , at 100 K from|Gurvich et al. (19891 



we find 7 = 1.33, 1.39, 1.33 and 1.30 for CH 4 , C 2 H 2 , C 2 H 4 and C 2 H 6 , respec- 
tively. Hydrodynamic atmospheric escape is quite sensitive to the value of 7. 
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The freezing out of the vibrational degrees of freedom, which increases the value 
of the adiabatic index, plays an important part in the theory's ability to match 
the spectroscopic observations from KBO surfaces. 
[Figure 1] 

In fig.l we plot the bulk density of four Kuiper belt objects as a function of 
their diameter. The error ellipses around each object indicate the uncertainties 
in the measured masses and radii of these objects. Also shown is the relation 
between bulk density and D max (from eq.[l]) for TJ, = 100 K for the four hy- 
drocarbons of interest. All objects to the left of a given curve will allow that 
particular gaseous species to flow hydrodynamically along an adiabat into space. 
All objects to the right of a given curve have a strong enough gravitational field 
to confine that particular gas which, upon re-cooling, ought to appear in the 
object's surface spectrum. The curves for C2H6 and C2H4 fall one on top of the 
other and are almost indistinguishable. We particularly chose intermediate size 
KBOs whose size and mass density are known, to a varying degree of accuracy, 
and can be compared with our theoretical predictions. In addition, the objects 
selected all show traces of crystalline water ice. Haumea presents a special case 
because of its rapid rotation and non-spherical shape, and we discuss it in more 
detail below. The error ellipse in this case was calculated using the diameter of 
a sphere with the same volume, with the uncertainties derived from conditions 



on an oblate spheroid as given by (Rabinowitz et al. 2006). 
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4.1 Slowly Rotating Bodies 



Orcus - The diameter and mass are taken from (Brown et al. 2010) and (Lim 



et al. 2010). According to (Delsanti et al., 2010) crystalline water ice is unam 



biguously detected. As pointed out by these authors, the noise in the data for 
Orcus is high and Hapke models do not give a unique solution. Therefore one 
cannot firmly exclude the presence of ammonia and methane on the surface. 
The observed spectrum is best fit by a mixture of water ice (mostly crystalline), 
ammonium (NH+) and ethane. C 2 H 2 and C 2 H 4 were not searched for since 
ethane was considered more likely, being an irradiation product of methane. It 
can be seen from fig. [T] that, according to our model, Orcus is not expected to 
retain C 2 Hg and methane on its surface. Higher molecular weight hydrocarbons 
can be present, however, and might explain the weak bands observed in Orcus' 
spectrum. 

Quaoar - Besides crystalline water ice, the presence of methane ice is also 
confirmed, where addition of ethane improves the quality of the fit to the spec- 



trum [see ( Guilbert et al. 2009 ) and references cited therein] . This is in complete 



agreement with our model, as shown in fig. 1. Quaoar's diameter is 890 ± 70 km 
and with the help of its satellite, Weywot, its density was determined to be 
4.2 ± 1.3 g cm -3 (Fraser and Brown 2010). According to eq. 16, Quaoar must 



have a mass density of/9&>3.8±0.6g cm -3 in order to show traces of methane 
ice on its surface. The uncertainty in this result stems from the uncertainty 
in Quaoar's reported diameter. We see that the exceptionally high density of 
this object could have been predicted without any measurement concerning its 
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satellite, Weywot. 

Charon - Occupation measurements give Charon's diameter as 1207±3 km 



and an average mass density of 1.71 ± 0.08 g cm" J (Sicardy et al. 2006). For 
these values, our model predicts no methane should be visible in the surface 
spectra, although C2H6, as well as C2H2 and C2H4 could be present. The 



surface spectrum of Charon was investigated by Cook et al. (2007) who found 
evidence for crystalline water ice. Adding methane ice to their model ice mixture 
did not improve the fit to the observed surface spectra of Charon. We found no 
mention that C2H6 was searched for. 

4.2 Rapid Rotators 

Haumea - Haumea presents a particularly interesting case. It is a fast rotat- 



ing elongated body that probably suffered a giant impact (Brown et al. 2007 



Delsanti et al. 2010). Pinilla-Alonso et al. (2009) showed that its surface is best 



described by amorphous and crystalline water ice and that it is probably de- 
pleted of volatile ices. As Haumea is a rapid rotator our model does not strictly 
apply, since we have assumed that the effect of the centrifugal acceleration im- 
parted by the rotating body on the escaping gas is negligible. For Haumea the 
rapid rotation effectively lowers the depth of the gravitational well. 

Since a rapidly rotating body will deviate from spherical symmetry and there 
is likely to be a superposition of tangential flow on a rapid radial flow, we need 



to modify the theory of coronal expansion as developed by E.N. Parker (Parker 



1963) for non-radial streamlines by adding the effect of centrifugal acceleration. 
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The flow velocity, v, along the streamline obeys: 

dv 1 dp 9$ 
V ~dl + ~ P ~dl + ~dl + ~dT = ° 

where I is the length along the streamline, p and p are the gas mass density and 
pressure respectively, $ and \& are the gravitational and centrifugal potentials 
external to the body. Assuming a polytropic relation between the gas pressure 
and density, with a polytropic index 7, one obtains: 

d ( l 2 7 Po -v-i 



pf-i + $ + * = (3) 



9/ V2^ ' 7-lp3' 
The subscript zero denotes the surface value of the given parameter. The ex- 
pression in parenthesis will be constant along the streamline. 



As noted by Parker (1963), conservation of matter along the streamline may 
be expressed as: 

p{l)v{l)A{l) = PqVqAq 

where A is the cross section of the flow tube. Combining this with eq.[3] and 
defining the relative cross section F = A/Aq gives: 

\ad + ^ r^r 1 + *w + = \ v i h — ^— — + * + * (4) 

2 7 — I po \vr J 2 7 — 1 po 

This is Parker's coronal flow equation including the centrifugal potential. 

If the streamlines are indeed non-radial then I, the distance along the stream- 
line, cannot be approximated by the radial coordinate r. A more natural choice 
for the independent spatial coordinate is the gravitational potential itself. Fol- 
lowing Parker, and including some changes due to the introduction of the cen- 
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trifugal acceleration, we define the following dimensionless parameters: 

v _j_ 

£ = — X'< +1 
vo 

= 7-1 PqVq 
27 p Q 

$ 2 

<j> = -2^x^ +1 

* 2 

With the help of these dimensionless variables eq.|3]can be rewritten as 

e 2 (<f>) + [F^)^)] 1 " 7 -<$>- i>(<p) = x^t + x ^ - O - fa (5) 

Let us examine the asymptotic behavior of this last equation along the 
streamline as (f> — > 0. It is reasonable to assume that far from the body stream- 
lines diverge, causing F to diverge. In addition, ip will vanish. After some 



algebra (see Levi and Podolak 20091 we again end up with two branches for 
the asymptotic solution far from the body, an upper branch and a lower branch, 
given to zero order by: 

Supper ~ \J X^ 1 + - 00 - 

eiower~{x^+X^-(t>o-^o) 1 "'F^ 1 (6) 
These are the same as those given by Parker, but with the addition of the fa 
term. Several complexities are introduced due to the streamlines' non-radial 
form that have to do with the way the flow transforms from subsonic to super- 
sonic flow. Nonetheless it is clear from the asymptotes that the establishment 
of a steady-state solution far from the body requires: 

x 4t + x ^i -<f> Q -fa > (7) 
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Reverting back to our original notation this may be written as: 



2(1-7) 
J_ Q„, 7+1 



+ 2v 



1 Po 



> 



(8) 



.7 - 1 Po 

In the limiting case where the body is massive enough to prevent the streamline 
from forming (i.e. vq — > 0) the inequality turns into an equality and since 7 > 1 
this reduces to the condition: 



7 kT Q 



*o + *o = 



(9) 



7 - 1 m g 

In cases where the centrifugal acceleration is of importance the object ro- 
tates quickly enough so that it deviates (perhaps considerably) from spherical 
symmetry. Let us assume the body behaves as a liquid rotator in hydrostatic 
equilibrium resulting in a problem independent of longitude. The shape of an 



ellipsoid of revolution may be formulated as (see, e.g. Zharkov and Trubitsyn 



1978) 



cos 2 9 



sin 



(10) 



.(1-e) 2 

Here ro is the radial coordinate at the body's surface, 9 is the colatitude, a is 
the equatorial radius, and e is the flattening parameter that, given the polar 
radius 6, is defined by 

a — b 



The gravitational potential on the surface of the ellipsoid may be expanded 
using gravitational moments as follows: 

GM 



$ = - 



ro(0) 



(11) 
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where G is the universal gravitational constant, M is the body's mass, J n are 
the gravitational moments and P n are Legendre polynomials. The centrifugal 
potential on the body's surface will be 



* = ~ny sm 2 9 



(12) 



where f^o is the object's angular velocity. Inserting all of this into eq. [9] results 
in the following restriction on the flow: 



7 kT GM / cos 2 



7- 1 m g 



(1-e) 2 
1 



cos 9 



n— 1 N v ' 



V&a 2 sin 2 I 



cos 2 9 



sin 2 6> 



sin 2 J 2 „P 2 (cos 0) 



(13) 



In the case where the body does not rotate eq. [13] reduces to the restriction given 



Levi and Podolak (2009) 



We now consider the case of Haumea and a possible methane atmosphere. It 
is not clear whether Haumea is an elongated ellipsoid or a flattened spheroid (see 



Rabinowitz et al. 2006). Ragozzine and Brown (2009) have adopted a tri-axial 



model. Nonetheless, we shall assume that the body can be approximated by a 
Maclaurin spheroid and apply eq.[l3] While this may not be strictly correct, it 
is useful in exploring some surprising consequences of rapid rotation for atmo- 



spheric escape. The body parameters are taken from Rabinowitz et al. (2006). 
We take a mass of M = 4.21 x 10 24 g, an average mass density of p w 2.6 g cm -3 , 
a polar radius of ~ 500 km, and an equatorial radius of ~ 800 km. The rotation 
period is taken to be 3.9154 hr. The quadrupole moment, J 2 , is not known. 



Ragozzine and Brown (2009) argued that J 2 ought be large due to Haumea's 



high amplitude light curve and adopted a value of 0.244 for their tri-axial model. 
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In order to understand the co-latitudinal restriction on the flow we examine 
a simple scenario; that of a constant density body. We define a rotational 
parameter, q, as: 



m 2 

4nGp 



0.2723 



(14) 



According to Zharkov and Trubitsyn (1978) we have for the constant density 
model: 



1 1 + ^ 



925 



15 
b 



1 



0.3801 



1 - e 

5 25 



806 \km] 



J 2 =2<7|l- n <T 



15 



98 



28 



125 



0.1255 



-0.032 



3 ~ 0.015 



(15) 



In fig. |] we plot the LHS of eq.[L3j the escape criterion, from the north pole 
(8 = 0) to the equator [6 = 7r/2) for a body with uniform density. Adiabatic 
flow is allowed only for colatitudes for which the escape criterion is positive. 
Thus, for the conditions we have assumed, adiabatic flow cannot be supported 
at a surface temperature of 130 K. However, if the surface could be heated to 
160 K, adiabatic streamlines may form for 9 > 0.427T, corresponding to latitudes 
in the range [—15° 15°]. Surface temperatures higher than about 220 K are 
required to support a global adiabatic escape. 

Taking the body's mass and volume as constant while varying J% has the 
effect of changing the body's internal density distribution. The effect of such 
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a variation on the escape criterion is shown in fig. [3j assuming a surface tem- 
perature of 150 K. The solid curves show the value of the escape criterion for 
three cases: J2 = 0.05, 0.1, and 0.2. Since the object's mass is derived inde- 
pendently of its shape we also consider a 10% increase in the equatorial radius 
(from 800 km to 880 km) while keeping the same mass. This is shown for the 
same three choices of J2 by the dashed curves. For all the cases presented in the 
figure adiabatic streamlines cannot be supported in mid-latitudes. For values of 
the quadrupole moment lower than the value corresponding to uniform density 
(~ 0.1), adiabatic streamlines can be supported only in a region around the 
equator, for the case of J2 =0.1 the latitudinal range allowing for adiabatic gas 
escape will increase from [—10° 10°] to [—20° 20°] when the equatorial radius 
is increased by 10%. This is due to the increase in centrifugal acceleration. 

For J2 = 0.2 adiabatic streamlines can be supported from the polar region, 
while only for the case of the larger equatorial radius there will be equatorial 
flow as well. This surprising behavior may be explained as follows: A low 
enough value of J2 represents a differentiated body with a dense core. Since 
the equatorial dimension is larger than the polar dimension, an observer at the 
pole will be closer to the core, and feel a stronger attraction from it, while 
the gravitational attraction from the lobes of the oblate mantle will cancel to 
some extent due to symmetry. Therefore, a denser core will result in a higher 
gravitational attraction at the poles and gas escape will be confined to the 
equator. Increasing the quadrupole moment means that more of the dense 
material is in the object's envelope rather than the core. This decreases the 
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gravitational attraction at the pole while somewhat increasing the gravitational 
attraction at the equator. As a result polar atmospheric escape is easier. In 
contrast to the poles, no cancellation due to symmetry occurs at mid-latitudes 
so the gravitational attraction there will remain high. This can be seen in 
fig. [4j where we have plotted the total (gravitational + centrifugal) potential 
as a function of co-latitude for the same three choices of J2. The behavior we 
describe is clearly seen. 

As methane sublimates from a given area it will flow (restricted by the escape 
criterion) along a streamline into space, clearing that region of methane ice. In 
fig. [5] we present the surface fraction evacuated as a function of surface temper- 
ature, for the three gravitational quadrupole moments and two equatorial radii 
considered in fig. [3] Due the object's geometry, a given latitudinal range near 
the equator corresponds to a larger surface fraction than would the same range 
near the pole. For the case of J2 = 0.05 the flow begins near the equator at low 
temperature and approaches polar regions monotonically for higher surface tem- 
peratures. This results in a rapid increase in surface fraction evacuated. Similar 
behavior is seen for J2 = 0.1 up to a certain high surface temperature where 
polar escape begins as well making the surface fraction evacuated increase even 
more rapidly. For J 2 = 0.2 flow begins from polar regions resulting in a rather 
slow increase of surface fraction evacuated per increase of surface temperature 
up to the temperature for which equatorial escape begins to take effect. 

From fig. [5] it is seen that at surface temperatures about half to two thirds 
of those required for global atmospheric escape, streamlines evacuating subli- 
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mating gases may already emanate from a large fraction (tens of percent) of 
the body's surface. If the axis of Haumea's rotation is perpendicular to our 



line of sight as suggested by Ragozzine and Brown ( 2009 ) and Rabinowitz et al. 



( 2006 ) then the equatorial region carries more weight observationally. In such 
a scenario the value we estimated for the surface fraction evacuated will be a 
lower limit to the value observed from Earth. For a = 800 km a high depletion 
of surface methane would suggest a differentiated body since surface temper- 
atures, even during extreme heating, do not go above 150 K. For the case of 
a = 880 km the same would be true if the surface temperatures do not go above 
130 K. 



5 Discussion 

As we have shown, the surface spectra of the four intermediate-size KBOs dis- 
cussed can be understood in the context of surface heating followed by rapid, 
hydrodynamic escape. The major advantage of this model is that very simple 
relations can be obtained between detection of certain volatiles and the body's 
diameter and density. We suggest the following scheme for bodies that show 
evidence of crystalline water ice. 

Given a body with a measured diameter D ± AD, that shows traces of CH4 
ice on its surface, its average mass density will obey: 



p b > 2.98 
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Since the curves for C2H6, C2H2 and C2H4 are almost indistinguishable, as seen 
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in fig. 1, we suggest that, if either one of these three hydrocarbons is observed 
while no methane ice is observed then: 
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(17) 



In case none of the four hydrocarbons (CH4, C2H2, C2H4 and C2H6) is observed 
we expect: 



Pb < 1.71 
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(18) 



As an interesting test case, let us consider the object 2003 AZ§4. Crystalline 



water ice was detected on its surface (Guilbert et al. 2009). Some discrepancy 



exists concerning its diameter. Stansberry et al. ( 2008 ) give a value of 686 ± 



95 km whereas Miiller et al. ( 2010 ) report a value of 896 ± 55 km. If methane ice 
will be discovered on its surface then, according to our model, its density would 
have to be higher than 6.3 ± 1.7 or 3.7 ± 0.5 g cm -3 , depending on its diameter. 
If CH4 ice is not detected but either C2H2, C2H4, or C2H6 are observed, the 
minimum density must be either 3.6 ± 1.0 or 2.1 ± 0.3 g cm -3 , respectively. If 
none of the four hydrocarbons is seen on the surface then the latter values are 
upper limits to the density. If the larger (smaller) estimate for the diameter is 
confirmed and methane (ethane) ice is detected then our model predicts this 
body has a Quaoar-like mass density, implying another KBO is asteroid-like, 
covered by a thin sheet of ice. This would have important implications for the 
way we view KBOs. If the body is rapidly rotating and has a large oblateness, 
as is the case with Haumea, observations of the surface coverage of volatile ice 
can be used to estimate the amount of differentiation the body has undergone. 
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Figure Captions 

Figure 1 - Atmospheric escape regimes for four hydrocarbon species compared 
with four KBOs. Error ellipses indicating the most likely uncertainties in the 
observed parameters are given. Being to the left (right) of a given curve means 
the gaseous species can (cannot) escape hydrodynamically. Haumea is a special 
case because of its non-spherical shape and rapid rotation, and is discussed sep- 
arately in the text. 



Figure 2 - Escape criterion (LHS of eq. 13 1 as a function of co-latitude for differ 



ent values of the surface temperature. Hydrodynamic escape occurs only when 
the escape criterion is positive (see text). 

Figure 3 - Escape criterion (LHS of eq.[l3| as a function of co-latitude for dif- 
ferent values of J2 and a surface temperature of 150 K. Solid curves are for an 
assumed equatorial radius of 800 km and dashed curves are for an assumed ra- 
dius of 880 km. 

Figure 4 - Total potential (gravitational plus centrifugal) on the surface of an 
ellipsoid of revolution as a function of co-latitude for different values of J2 . 
Figure 5 - Fraction of the body's surface from which adiabatic streamlines can 
be energetically supported as a function of surface temperature for different as- 
sumed values of J^. The solid curves are for an assumed equatorial radius of 
800 km. The dashed curves are for a radius of 880 km. 
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Flow Regimes for Various Hydrocarbons 
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Figure 1: Levi and Podolak, Estimating KBO densities 
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Colatitude extent of escape (uniform density model) 
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Figure 2: Levi and Podolak, Estimating KBO densities 



34 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

9[pi] 



Figure 3: Levi and Podolak, Estimating KBO densities 
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Figure 4: Levi and Podolak, Estimating KBO densities 
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Figure 5: Levi and Podolak, Estimating KBO densities 
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